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' Abstract 

New solutions on the scalar - isoscalar tttt phase shifts are analysed 
^ ■ together with previous KK results using a separable potential model of 

^ . three coupled channels (vrvr , KK and an effective 27r27r system). Model 

parameters are fitted to two sets of solutions obtained in a recent analysis of 
the CERN-Cracow-Munich measurements of the iT~pi — > iT~^7r~n reaction 
on a polarized target. A relatively narrow (90 - 180 MeV) scalar resonance 
/o(1400 — 1460) is found, in contrast to a much broader (F ^ 500 MeV) 
state emerging from the analysis of previous unpolarized target data. 

The spectrum of scalar mesons is still not well known [|, |]. The J^^ = 0++ 
meson nonet is not free from ambiguities and the internal structure of scalar 
mesons is rather controversial. Recently, many papers have been devoted to 
the study of the lowest scalar glueball [3-6]. In these articles properties of the 
/o(1500) resonance, observed recently by the Crystal Barrel Collaboration in pp 
annihilations JT], are analysed and interpreted. Furthermore, lattice QCD 
calculations predict the mass of the lowest scalar glueball to fall between 1500 and 
1700 MeV IP, |10|- One should also mention increasing evidence for a broad scalar- 
isoscalar resonance a [11-14] below 1000 MeV. This state reappeared as /o(400 — 
1200) in the last edition of the Particle Data Group [|l|], after 22 years of absence. 
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Well established resonance /o(980) lies at the KK threshold. Its interpretation, 



however, is a subject of lively discusions, as can be seen for example in Ref. [jT5| 
and articles cited therein. 

In Ref. |jTl| the scalar meson spectrum was studied in terms of a relativistic 
TTTT and KK coupled channel model from the tttt threshold up to 1400 MeV. 
The phenomenological parameters were constrained by fitting the 5'-wave data 
extracted from the experimental cross sections on the tt+tt^ production by tt" 
scattering on unpolarized hydrogen target Further constraints were imposed 
by the KK phase shift analysis of Ref. 



Recently, authors of Ref. [|l8l have analysed data obtained on a polarized 



target by the CERN-Cracow-Munich group for the 7r~p| — 7r"'"7r~n reaction [|1^ 
Separation of the vr and Oi exchange amplitudes in this reaction was then possible 
for the first time, using assumptions much weaker than in all previous analyses. 
From a set of four solutions for the isoscalar S'-wave phase shifts up to 1600 
MeV, two of them ("down-fiat" and "up-flat") satisfy the unitarity constraint. 
The "down-fiat" solution is in good agreement with the former solution of Ref. 
TB[] up to 1400 MeV. Above 1400 MeV one observes an increase of the phase shifts 



and larger inelasticies than those in Ref. |T^. This could be a manifestation of 
the presence of scalar mesons /o(1370) or /o(1500) in that energy range. There, a 
strong four-pion production has been observed in different experiments [7,8,20- 
22] . This provides a compelling argument to take into account the 47r channel. In 
this channel there is some evidence of clustering into a a or pp pairs [6-9] . In [^ 



an attempt was made to describe the tttt S'-wave in the 1000 to 2000 MeV region 



by generalization of the non-relativistic 2-channel model ||2^ to three coupled 
channels. 

In the present paper we extend the isospin S'-wave relativistic 2-channel 
model of Ref. |Tl| by adding to its tttt and KK channels an effective third 
coupled channel, here called aa. 

We consider the S'-wave scattering and transition reactions between three 
coupled channels of meson pairs labelled 1, 2 and 3. Reaction amplitudes T 
satisfy a system of coupled channel Lippmann-Schwinger equations ||11||: 



< p|T|q >=< p|y|q > + / ,f < p|V"|u >< 'a\G\\i >< u|T|q >, (1) 

where V, G and T are 3x3 matrices, V is the interaction matrix and G is the 
diagonal matrix of channel propagators: 



In Eq. (0) E is the total energy, Ej = yfn'j + and rrij is the meson mass 
in channel j. We consider meson pairs of same mass and momentum u in their 
centre of mass system. 
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We choose a separable form of the interaction: 



< P|K7,i|q >= X^-^aT.i fi'a,i(p) fi'7,j(q)' ",7 = 1,2,3, (3) 
where A^^ are couphng constants and 



[Att 



are form factors which depend on the relative centre of mass meson momenta p in 
the final channel or q in the initial channel. In the vrvr channel (j = 1) we choose 
a rank-2 separable potential (n = 2) and in the other channels, i.e. KK {j = 2) 
and aa {j = 3), a rank-1 potential (ra = 1). 

Following the formalism developed in Ref. we can solve the system of 
equations (|l]) - (|^), which leads to the Jost function 

Z}(ki,k2,k3) = rfet(l-A/), (5) 

where A and / are symmetric 4x4 matrices 
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(6) 



and 



(7) 



In (^) and (|^ k^ are the on-shell momenta connected with the total energy by 
E = 2Jkl + ml = 2Jkl + ml = 2^1 + ml (8 



Analytical expressions for la^ijika) integrals are given in Appendix A of Ref. [|T^ . 
Altogether, this model has 13 parameters: 9 coupling constants Aq,^ j and 4 range 
parameters I3a,j (for a = 1 j = 1 or 2 and for a ^ 1 j = 1). In the calculations 
we shall use dimensionless coupling constants defined as 



A, 



A 



"7: J 



"7. J 



2 iPajPyj) 



(9) 



From now on, we will omit index j for Xayj if both a and 7 are different from 1 
and for Pa,j if a 7^ 1. 
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S'-matrix elements Sap = 1,2,3) can be written in terms of the Jost 

function of different arguments, for example 



^11 



D(-ki,k2,k3 
/^(ki,k2,k3) 



(10) 
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Further expressions for some other matrix elements can be found in Ref. 
The model satisfies the unitarity condition S^S = 1. Diagonal matrix elements 
are parametrized as 



J = 1,2,3, (11) 

where r]j and 6j are channel j inelasticities and phase shifts, respectively. Expres- 
sions for nondiagonal elements can be found in ||25|. Some of the S-matrix poles 
in the complex energy plane can be interpreted as resonances. They correspond 
to the zeroes of the Jost function -D(ki, k2, ka) as can be seen from (plOl). 

We fit the existing experimental results on the tttt S'-wave isoscalar phase 
shifts together with inelasticity in the vrvr channel and with the KK phase shifts. 
In this analysis we extend the 2-channel model developed in Ref. |jll| since new 
data have been recently obtained |]TB|] in a wider mass range (up to 1600 MeV). 



We choose a method based on the x fit: 



2 2,2,2 
X =Xn + XnK + X^ 



rp 



(12) 



XnK 



Xn 



E 

1=1 



Xrj 



E 



2A(5^^ + 6kk) 



exp \ ■ 
iinr_ I 




(13) 



(14) 



(15) 



where 5,r7r,5x_ft: are phase shifts in channels 1 and 2, A(TT,,AaK are the corre- 
sponding experimental errors and r].jrn , Arj^m are inelasticities and inelasticity er- 
rors, respectively. Superscripts "th" or "exp" refer to theoretical or experimental 
values. In Eqs. (|1^) to (plSf ) Nt,,N.„k and A^,, are the numbers of experimental 
points. 

Below 600 MeV we use data from the KeA decay |2^ and from Refs. 
and ||28|. Above 600 MeV we use the "down- flat" and "up-flat" solutions of the 
analysis of Ref. |T3 on a polarized target. Constraints in the K^K~ channel 



are also needed. Therefore, we have used the results of the analysis of reactions 
Ti^p — s> K^K"n and vr^n K^K^p |]T^, although targets were unpolarized 
there. 

We begin with a simpler ^^-minimization by considering only the vrvr and KK 
2-channel case. Starting parameters of the interaction potentials have been taken 



from 2-channel fits obtained in ||29l]. In this case there are 8 free parameters, as 
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indicated in Table |l]. The corresponding fits for 5^,^, "/^Trvr and ip^^K = ^mr + ^kk, 
plotted as dotted lines, are compared to the experiment in Figs. 1, 2 and 3. Above 
1400 MeV both "down-fiat" and "up-fiat" data indicate a decrease towards small 
values of (?7 ~ 0.6 to 0.7), albeit with large errors. Furthermore, it can be 
seen in Fig. 2 that the errors of r] are much larger than dispersion of their values. 
Therefore, in order to reproduce more easily this trend of 77 as a function of energy, 
for most of the data points we reduce in Eq. (15) the experimental rj errors to 
0.1. For the "down-fiat" solution we keep, however, the original errors at 1310, 
1470 and 1570 MeV since inelasticities at these energies are above or close to 1. 
For the last point at 1590 MeV we increase the error from 0.41 to 1 since rj is as 
large as 1.52. Similarly, for the "up-fiat" solution we do not change the errors in 
the range between 1350 and 1450 MeV, however, at 1310, 1470 and 1590 MeV we 
increase the errors to 0.8 since rj values at those energies are substantially higher 
than 1. 

We found that in the 2-channel case for the "down-fiat" solution it was not 
possible to obtain the rj values with such modified errors between 1400 to 1600 
MeV (this is slightly different for the "up-flat" solution as can be seen in Fig. 
2b). In the 3-channel model, however, we can get a substantial decrease of 77 
above 1400 MeV (see Fig. 2a). In order to achieve this behaviour, couplings 
between the vrvr and aa or KK channels should be sizable. The corresponding 
14 free parameters including the mass for two of our best flts are given in 
Table |I|, and their in Table |^. The 3-channel flts for the " down-flat" data are 
labelled by A and B, and those corresponding to the "up-flat" data - by C and 
D, respectively. The flts favour 7713 masses in the range of 675 to 733 MeV. The 
values of S.,riT,VTnT and ipj^K are drawn as thick and thin solid lines in Figs. 1-3 
and compared to data. In Fig. 1 we have shown flts A and C only since energy 
dependence of flt B is very close to that of A, and that of D - very close to C. 
In Fig. 1 one can already notice a better agreement with data of the 3-channel 
model in comparison with the 2-channel one. This is especially well visible in Fig. 
lb. The main difference between the 2- and 3-channel flts lies in r] above 1400 
MeV, where the opening of the aa channel leads to a fast decrease of inelasticity 
parameters (see Fig. 2a). Let us also note an improvement in ip^^K over the 1000 
to 1200 MeV range, as can be seen in Fig. 3 and in Table ^ with a better xix- 

Fits of similarly good quality were obtained with very different physical pa- 
rameters in the KK and aa channels. For example, in the 2-channel flts and 
in flt A, the KK interaction is attractive while in the other cases it is repulsive 
(Table |l|). Similarly, interchannel couplings are very different in both cases. In 
flt A we see rather strong, and in flt C very strong tttt to aa and KK to aa cou- 
plings, while in the B and D cases the vrvr - KK coupling, Ai2,2, is particularly 
strong. Lack of a sufficient number of observables and/or experimental precision, 
in particular in the effective 47r channel, leads to the existence of several good sets 
of model parameters with quite different channel and interchannel interactions. 
Other 3-channel flts of reasonable quality with less than 14 parameters can be 
also obtained. 
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We have studied positions of the S'-matrix poles in the complex energy plane 
[E = M—iV/2). For the 3-channel model there are 8 different sheets which corre- 
spond to different signs of imaginary parts of channel momenta {Impi, /mp2, I^P?) 

For example, on the sheet denoted by ( ) all imaginary parts are negative. In 

the case of the 2-channel model there are only 4 sheets labelled by signs of Impi 
and Imp2. Resonance parameters predicted by the 2- and 3-channel models are 
summarized in Table 3. At low energy we find a very broad /o(500) resonance 
(also called cr meson). Since the "up-flat" data indicate a faster increase of the 
TTTT phase shifts with energy than the "down-fiat" data, the a meson appears 
in the first case at a mass higher by about 40 MeV and with a width smaller 
by about 50 MeV. In the three channel fits the /o(980) resonance is seen in the 
vicinity of the KK threshold with a width of about 60 to 70 MeV. The relatively 
narrow state /o(1400) appears in the 3-channel fits (see Table Its mass varies 
from about 1400 MeV to 1460 MeV. For the "down-flat" fits this resonance is 

narrower (F ^ 100 MeV) on sheet ( h) than on sheet ( ). This width is 

very close to values found by the Crystal Barrel Group (Refs. 0, 0) but the res- 
onance masses are smaller than their values close to about 1500 MeV. We should 
point out that a narrow (F = 65 ± 10 MeV) scalar resonance at M = 1445 ± 5 
MeV has been found by the WA91 group at CERN in central production of An 



in high-energy pp collisions |22 



Finally, let us mention an important difference between the fits presented 
here for data taken on a polarized target in comparison with those performed 
with data [0 obtained on a nonpolarized target. If we fit the data of [|T6| using 



our 3-channel model, we obtain a very wide resonance at M = 1521 MeV of 
width 503 MeV. This means that the recent analysis of data [1^ supplies some 
new information on the /o(1400) meson. 

In conclusion, we have analysed the isoscalar S-wave tttt and KK scattering 



using the CERN-Cracow-Munich data |jT8| together with the data of [17,26-28 



in the framework of the 2- and the 3-channel models of meson-meson scattering. 
All fits of the phase shifts analysis indicate presence of a relatively narrow 
(90 - 180 MeV) scalar resonance of mass 1400 - 1460 MeV. This resonance is 
quite compatible with recent observations of a possible scalar glueball at 1500 
MeV with a width of 100 MeV. 

We are indebted to D. V. Bugg for useful correspondence and to J. Kisiel, V. 
E. Markushin and K. Rybicki for fruitful discussions. 
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Figure 1: Energy dependence of vrvr phase shifts: a) fit to "down-flat" data 
of [18], thick solid line corresponds to fit A of the 3-channel model and dotted 
line to the 2-channel model fit; b) fit to "up-flat" data of [18], thick solid line 
corresponds to fit C of the 3-channel model and dotted hne to the 2-channel 
model fit. 
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Figure 2: Energy dependence of inelasticity parameter 7]^^,^: a) fit to "down- 
flat" data of [18], thick solid line corresponds to fit A, thin solid line to fit B and 
dotted line to the 2-channel model fit; b) fit to "up-fiat" data of [18], thick solid 
fine corresponds to fit C, thin solid fine to fit D and dotted fine to the 2-channel 
model fit. 
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Figure 3: Energy dependence of phase shifts sum ip^^K = ^tt-k + ^kk- ^) 
"down-fiat" data of [18], thiclc solid line corresponds to fit A, thin solid line to fit 
B and dotted line to the 2-channel model fit; b) fit to "up-fiat" data, thick solid 
fine corresponds to fit C, thin solid fine to fit D and dotted fine to the 2-channel 
model fit. Data for ip^^K are taken from [17]. 
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Table 1: Separable interaction parameters for 2- and 3-channel model fits to the 
"down-flat" and "up-flat" data from [18]. Values of (3 and 7713 are given in GeV. 



Data 


down-flat 


monpl 




3-channel 


fit 

lit 




A 


B 


Aii,i 


-.14258 X 10-^ 


-.29975 X 10-2 


-.52138 X 10-2 


All,2 


-.18895 


-.10844 


-.10552 


A22 


-.49106 


-.39304 


3.1637 


A33 





-.17447 X 10-2 


-.45719 X 10-1 


Al2,l 


.27736 X 10-5 


.12039 X 10-2 


.10685 X 10-1 


Al2,2 


.43637 X 10-1 


-.11333 


-.80626 


Al3,l 





-.25841 X 10-2 


.14544 X 10-* 


Al3,2 





.39924 


.21878 X 10-2 


A23 





-.57955 


-.17515 X 10-1 


A,i 


3.0233 X 10^ 


1.426 X 102 


.81615 X 102 


/9l,2 


1.0922 


.92335 


.85776 


(^2 


2.2941 


1.4959 


.47403 


P3 




1.3676 


.45357 X 102 


772,3 




.70 


.67510 


Data 


up-flat 


mnn p1 


2-cliannel 


3-channel 


fit 




C 


D 




-.25550 X 10-* 


-.24649 X 10-2 


-.13063 X 10-2 


All,2 


-.15677 


-.73339 X 10-1 


-.94483 X 10-1 


A22 


-.66905 


.19475 


3.6391 


A33 





-.45073 X 10-2 


-1.4588 


Al2,l 


-.23243 X 10-5 


.12710 X 10-2 


.10185 X 10-2 


Al2,2 


.74350 X 10-1 


-.32171 


-.85407 


Al3,l 





-.52567 X 10-2 


-.18043 X 10-=^ 


Al3,2 





1.6134 


.17803 


A23 





-4.0462 


-.50484 




1.6878 X 10* 


1.7379 X 102 


3.2899 X 102 


A,2 


1.5882 


.92153 


.90097 




1.4710 


.69032 


.30733 






.23546 


1.1419 


ma 




.73252 


.68096 
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Table 2: Best values of 2- and 3-channel model fits for the "down-flat" and 
"up-flat" data of ref. [18]. Numbers of data points are indicated in parentheses, 
values are obtained when fitting with reduced 77 errors. 



Data 


down-flat 


up-Hat 


Model 


2-channel 


3-channel 


2-channel 


3-channel 


fit 




A 


B 




C 


D 


xl (65) 
xUk (21) 

xl (30) 
xlt (116) 


66.4 
26.3 
9.4 

102.1 


63.0 
15.9 
13.2 

92.1 


61.2 
9.7 
12.9 

83.8 


115.6 

27.0 

10.9 

153.5 


88.0 
13.1 
14.1 

115.2 


85.3 
13.2 
14.6 

113.1 


X^'' (30) 
Xtot^ (116) 


100.6 
193.2 


36.7 
115.6 


29.3 
100.1 


43.3 
186.0 


36.2 
137.3 


37.7 
136.2 



Table 3: Masses and widths of resonances found for the 2- and 3-channel fits to 
the "down-flat" and "up-flat" data 



2-channel model 



pole 


sheet 


down-flat 


up-flat 






M (MeV) 


r (MeV) 


M (MeV) 


r (MeV) 


/o(500) (a) 


-+ 


524.0 


513.6 


592.9 


401.8 


/o(980) 


-+ 


993.4 


79.3 


1015.0 


101.4 


/o(1400) 




1434.6 


167.6 


1429.7 


179.2 


3-channel model 






down-flat 


pole 


sheet 


A 


B 






M (MoV) 


r (MeV) 


M (MeV) 


r (MeV) 


/o(500) (a) 


- + + 


518.1 


521.4 


511.8 


532.6 


/o(980) 


- + + 


989.0 


62.0 


992.4 


68.2 


/o(1400): 




1405.1 


147.8 


1411.5 


169.3 




h 


1456.4 


93.3 


1402.7 


108.2 








up- 


flat 




pole 


sheet 


C 


D 






M (MeV) 


r (MeV) 


M (MeV) 


r (MeV) 


/o(500) (a) 


- + + 


561.7 


466.8 


558.2 


477.6 


/o(980) 


- + + 


992.2 


67.6 


995.3 


71.7 


/o(1400): 




1407.0 


181.2 


1418.8 


179.3 




h 


1423.8 


177.4 


1416.5 


173.2 
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